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THE  CHARGE-EXCHANGE  NEUTRAL  PROBLEM  AND  THE 

CODE  NEUTRAL 


Introduction 

In  this  report  we  will  discuss  the  theoretical  aspects  and  the 

computer  modeling  of  the  charge-exchange  neutral  process.  Such  a 

mechanism  has  been  proposed  to  explain  the  premature  gap  closure  in 

pulsed  power  ion  diodes  as  well  as  in  the  Reflex  Switch  .  In  the 

reflex  switch  it  is  important  that  ve  understand  this  anamolous  gap 
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closure,  since  such  an  event  causes  reduced  power  to  the  load  .  In 
the  presentation  below  we  briefly  review  the  operation  of  the  reflex 

4 

switch  and  then  we  will  discuss  the  analytical  model  for  charge- 
exchange.  We  then  will  present  some  results  from  a  monte  carlo 
simulation  code. 


Reflex  Switch 

The  closed  mode  of  the  reflex  switch  can  be  seen  in  figure  1  .  An 

anode  foil,  denoted  here  in  the  text  by  A,  is  set  between  the  primary 

cathode  K1  and  the  electrically  floating  secondary  cathode  K2.  There 

is  also  an  applied  axial  magnetic  field  Bz,  which  we  will  not  concern 

ourselves  with  here.  In  the  closed  mode  of  the  switch,  electrons 

emmitted  from  K1  are  repelled  by  the  potential  in  the  A-K2  gap, 

causing  the  electrons  to  repeatedly  traverse  (reflex)  through  the 

anode  foil.  From  the  resulting  anode  plasma  an  ion  current  is  drawn, 

as  shown  in  the  figure.  As  these  ions  are  accelerated  through  the 

partially  ionized  anode  plasma,  they  can  charge-exchange,  creating  a 

flux  of  energetic  neutrals  into  the  gap.  The  neutrals,  by  electrical 

breakdown,  then  cause  an  impedance  collapse  in  the  anode-cathode  gap*. 

This  is  the  model  of  neutral-charge  exchange  as  mentioned  in  the 

introduction.  The  anode-K2  gap  spacing  determines  when  the  switch 

opens,  and  it  is  supposedly  plasma  motion  that  causes  this  gap  to 

short  and  stops  the  electrons  from  reflexing^.  It  is  also  possible 

that  neutrals  can  be  filling  this  gap  as  well,  making  the  switch  open 
o 

prematurely  .  We  will  now  go  into  more  detail  on  the  neutral  charge- 
exchange  problem  as  given  below. 

Manuscript  approved  August  31.  IVtSSi. 


1 


Neutral  Charge-Exchange  Hodel 

Prono  et  al.^  were  the  first  to  suggest  that  premature  gap  closure  was 

caused  by  neutrals  that  were  produced  from  ions  charge-exchanging  in  a 

background  medium.  The  physical  model  they  used  is  given  in  figure  2. 

There  is  an  ion  flux  at  x  =  0  into  a  partially  ionized  plasma,  which 

contains  mostly  neutrals  at  a  uniform  density  nA*  The  ions  charge- 

exchange  as  they  are  accelerated  through  the  neutral  gas,  by  the 

constant  electric  field  E.  In  the  microscopic  picture  of  the  charge- 

exchange  process,  an  ion,  given  by  the  "+",  picks  up  an  electron  from 

a  similiar  type  atom(s),  producing  an  energetic  neutral  n  with  the 

same  velocity  as  the  incident  ion  and  an  ion  with  zero  speed.  We  note 

that  this  newly  born  ion  will  then  gain  energy  from  the  field  E,  and 

can  charge-exchange  as  well.  Thus  one  ion  can  produce  many  neutrals. 

The  flux  of  neutrals  into  the  vacuum  gap,  given  by  G(x=x  ),  is  the 

max 

cause  of  the  gap  closure.  The  physical  model  of  the  charge-excahnge 
process,  as  given  here,  can  be  described  analytically  and  is  presented 
below. 


To  find  the  neutral  flux,  we  follow  the  analysis  given  by  Prono 
et  al.*.  Assuming  the  problem  described  above  is  in  a  steady  state  we 
can  solve  for  the  ion  distribution  function  f(x,v)  and  the  neutral 
distribution  function  g(x,v)  by  starting  with  the  one  dimensional 
Boltzmann  equation  given  by 

*  2s  If  ■  -v£  *  «v-*>£v<v-vf'  (i) 

v||  .  nA»vf  (2) 

where, 

f(v,x)dxdv  =  number  of  ions  (per  unit  area)  in  the  phase 
volume  dxdv 


g(x,v)dxdv  =  number 
volume 


of  neutrals 
dxdv 


(per  unit  area)  in  the  phase 
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q  =  charge  of  the  interacting  ion 
m  =  mass  of  the  ion 
E  =  electric  field 
nA  =  neutral  background  density 
v  =  velocity  in  x  direction 

a  =  charge  exchange  cross  section,  which  in  general 
depends  on  v 

By  a'  and  f'  in  equation  (1),  we  mean  f'  =  f(v' ,x)  and  a'  = 
c(v').  Also  in  (1)  we  have  assumed  that  v  is  always  greater  than 
zero  (no  reverse  flow),  and  initially,  we  assume  that  the  slow 
ion  produced  in  the  charge-exchange  process  has  a  very  small 
energy.  This  is  the  meaning  of  the  e  term  in  the  argument  of  the 
Dirac  delta  function.  So  the  first  term  on  the  right  hand  of 
equation  (1)  is  the  loss  of  ions  due  to  charge-exchange,  while 
the  second  term  denotes  ion  production.  While  in  equation  (2)  we 
see  that  the  neutrals  are  not  affected  by  the  electric  field,  and 
that  nA«rf  is  the  neutral  production  rate.  Also  note  that  in  the 
above  equations,  as  well  as  in  our  model  in  general,  we  have 
ignored  the  ionization  of  the  fast  neutrals  as  they  pass  thr  ugh 
the  background  medium.  Steady  state  conditions  imply  that,  the 
flux  of  ions  is  a  constant,  which  we  and  Prono  et  al.  denote  as 
F.  Thus  we  can  write,  assuming  the  ions  are  born  with  small 
velocity  e, 


F-I, 

dvvf (v,x) 

> 

(3) 

f(v,x=o)  =  | 

8(v-e) 

(4a) 

s 

< 

II 

o 

II 

o 

(4b) 

where  (3)  is  just  the  usual  definition  of  the  flux,  and  equation  (4a) 
is  the  assumed  ion  distribution  at  the  injection  point.  The  quantity 
e  is  used  to  avoid  mathematical  ambiguities  later  on  in  the  analysis 
and  will  eventually  be  set  to  zero1-  Equation  (4b)  states  that  there 
are  no  neutrals  created  at  the  injection  point.  We  can  Laplace 
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transform  equations  (1)  and  (2),  using  the  boundary  conditions  (4a) 
and  (4b),  to  find 


31 

3w  + 


o)I 


g  = 


6(w-e)mF 

qE 

n.  <rf 
A 

P 


n.  a'  f ' 
A 

mF 


(5) 

(6) 


with  the  usual  definition  of  the  Laplace  transform  being, 

r® 

I  =  f(v,p)  s  dx  e~px  f(v,x) 

“  o 


Note  that  in  equations  (5)  and  (6)  we  have  used  a  change  of  variable, 

2 

where  w  is  the  kinetic  energy,  w  =  mv  /2.  Solving  for  f,  after  some 
integrations  of  equation  (5),  see  Appendix  A,  we  find 


?(w)  =  f(w,p)  = 


f  H(V)  -£■ 

^  Oft 


V 

-  I  5e  (p  * 


,  -  I  5b  (p  4  V'W 

dw  e  Jo  * 


H(w) 


1  for  w  >  0 
0  for  w  <  0 


(7) 

where  H(w)  is  the  usual  Heaviside  step  function  and  we  have  let  e  go 
to  zero.  From  equations  (7)  and  (6)  we  see  that  once  we  have  the 
functional  form  of  the  charge-exchange  cross  section  a,  we  can  in 
principle  determine  the  ion  and  neutral  distribution  functions  by 
inverse  Laplace  transform.  However,  one  can  see  that  in  general, 
exact  solutions  will  be  very  difficult,  if  not  impossible,  to  obtain. 
At  this  point  we  diverge  from  the  analysis  of  Prono  et  al.,  and  we 
assume  that  o,  the  charge-exchange  cross  section,  is  constant  for  all 
w.  Even  though  this  assumption  is  less  physical  than  having  a  cross 
section  that  decreases  with  energy,  as  Prono  et  al.  did  ,  we  can  at 
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least  obtain  closed  analytical  expressions  for  f  and  g.  We  note  here 
that  if  the  average  ion  energy  is  much  less  than  the  "cut  off  energy" 
of  the  cross  section,  again  see  reference  (1),  then  a  constant  cross 
section  is  a  good  approximation.  Returning  to  the  problem  at  hand, 
using  equations  (6)  and  (7)  with  a  being  a  constant,  we  find  by 
inverse  Laplace  transform  that  the  distribution  functions  f  and  g  are 
given  by 


f(v,x) 


g(v,x) 


(8) 

(9) 


where  X  =  1/n^cr,  which  is  the  mean  free  path  for  charge-exchange.  As 
before,  H  is  the  Heaviside  step  function  as  defined  in  equation  (7), 
and  5  is  the  usual  Dirac  delta  function.  Taking  various  moments  of 
the  distribution  functions,  we  find 


vf(v,x)dv  = 
Jo 

J"£(»,x)^  .  F 

(10) 

rn<x>  - 

*  F-  (  l  ) 

(11) 

ao 

Fu(x)  '  JoW*<w,x»gS£.  . 

F' (qEX) (  1  -  e'(x/X)  j 

(12) 

00 

G„(x>  -  J  »r(w,x)^  . 

F- (qEX) (  (e~(x/X)-  1  )  +  (x/X)  ) 

(13) 

Equation  (10)  is  just  a  verification  of  equation  (3).  rn(x)  *s  the 
neutral  flux  profile,  which  for  our  case  happens  to  be  linear.  F  (x) 
and  Gy(x)  are  the  ion  and  neutral  energy  fluxes,  respectively.  Notice 
that  the  term,  qEX,  keeps  appearing  in  our  equations.  This  is  just 
the  average  energy  an  ion  gains  before  charge-exchanging.  Defining 
the  mean  particle  energy  as  Prono  et  al.  does''’,  we  find 
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p 

<  E.>  =  f-  =  (qEX) (  1  -  e_(x/X)  )  (14) 


G 

<E„>  '  -7  ' 

n 

(qEX) 

> 

)(e-(x/X) 

1 

+ 

1— » 

t _ l 

(15) 

both  cases,  as  x 

0 

t 

both 

<E . >  and 

1 

<En>  go  to 

zero,  as  one 

would  expect.  For  x  -»  °°,  both  expressions  go  to  the  value  qEX.  The 
analytical  solutions  (10)  through  (15)  are  very  useful  in  benchmarking 
a  computer  simulation  of  this  neutral  charge-exchange  model.  In  the 
next  section  we  will  compare  (9),  (10),  (11),  and  (15)  to  the  results 
from  a  computer  code  named  NEUTRAL. 


Computer  Simulation  and  Modeling 

The  code  NEUTRAL  was  developed  to  model  the  charge-exchange 
problem  as  mentioned  in  the  previous  section,  and  is  the  final  member 
in  a  series  of  codes  which  varied  in  their  degree  of  complexity.  The 
first  code  that  was  written  was  designed  to  simulate  ions  interacting 
with  a  background  medium  with  no  applied  electric  field  and  no  time 
dependence.  The  second  code  in  the  series  included  ion  charge¬ 
exchanging  in  an  applied  electric  field,  but  the  model  was  still  time 
independent.  The  final  member,  NEUTRAL,  includes  all  ot  the 
capabilities  of  the  previously  mentioned  codes,  and  is  time  dependent. 

NEUTRAL  uses  Monte  Carlo  techniques  to  model  the  charge-exchange 
process,  and  the  code  was  written  so  that  either  applied  or  self- 
consistent  electric  fields  can  be  used.  This  second  field  mode 
however,  has  not  been  tested  yet.  Our  main  concern  was  to  model  the 
problem  as  envisioned  by  Prono  et  al.^,  since  experimental 
measurements  of  relevant  quantities,  such  as  densities  and  electric 
fields,  are  practically  nonexistent.  The  code  works  by  imitating  the 
problem  as  pictured  in  figure  2.  Ions  are  injected  at  the  left  hand 
boundary  and  then  an  electric  field  is  determined.  Next  the  particles 
are  moved  in  position  and  velocity  space.  Then  the  code  decides 
wether  or  not  a  given  ion  undergoes  charge-exchange.  This  is  done  by 
randomly  generationg  a  number  between  0  and  1,  and  then  comparing  it 
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to  e-^k/X^  where  AL  is  the  distance  traveled  by  the  particle  in  a 
given  time  step.  As  before,  X  is  the  mean  free  path  for  charge- 
exchange.  If  the  random  number  is  greater  than  e  a  charge- 
exchange  event  event  is  said  to  have  occurred,  and  in  the  fast  ions 
place,  a  fast  neutral  with  the  same  velocity  is  created  and  tracked, 
while  another  ion  is  created  with  zero  speed.  One  can  see  that  this 
method  assumes  that  an  ion  interacts  only  once  in  a  given  time  step, 
which  is  valid  if  AL  <<  X.  This  method  is  also  well  suited  for  the 
case  where  X  depends  on  x  and  w,  through  nA  =  nA(x)  and  a  =  a(w).  So 
to  summarize  NEUTRAL,  for  every  time  step,  ions  are  injected,  electric 
fields  are  computed,  particles  are  pushed  in  x  and  v,  and  charge- 
exchange  events  are  randomly  generated.  We  now  will  present  some  of 
the  output  from  the  code  NEUTRAL  in  the  next  section. 


Simulation  Results 


As  mentioned  earlier  we  are  using  the  problem  posed  by  Prono  et 
al.1,  assuming  the  charge-exchange  cross  section  a  is  a  constant  in 
energy,  to  benchmark  the  code  NEUTRAL.  We  used  protons  as  our  ion 
source  and  had  hydrogen  as  our  background  gas.  Some  of  the  physical 
parameters  assumed  in  this  standard  (reference)  case  were 


thus, 

and 


n.  =  1.0  x  10  /cm 
A  o 

E  =  1  MV/cm  =  3.333333  x  10J  statvolts/cm 

q  =  4.8032  x  lO-1^*  statcoulombs 
-15  2 

a  =  2.0  x  10  ctn 
m  =  1.6726  x  10“24  g 


X  =  5  x  10'2 


cm 


NTS  = 


500 

400 


At  =  1.136  x  10 
to  =  5.0  x  10  ^ 


11 

cm 


sec 
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The  quantity  m  is  the  ion  mass,  which  for  this  case  is  just  the  proton 
mass,  and  q  is  the  charge  of  the  proton.  From  the  above  input 
parameters  it  would  not  be  surprising  if  we  said  that  NEUTRAL  does 
most  of  its  calculations  in  c.g.s.  units.  The  second  set  of  variables 
^inj’  ^t,  an<*  A*  are,  respectively,  the  number  of  ions  or 
particles  injected  per  time  step,  the  number  of  time  steps  run  in  the 
simulation,  the  time  step  in  seconds,  and  the  size  of  the  spatial 
mesh.  The  time  step  was  chosen  such  that  At  =  Ax/v  ,  where  v  is 
the  velocity  of  an  ion  if  it  traversed  the  entire  charge-exchange 
region  without  interacting.  The  quantity  xmax  is  the  position  of  the 
right  hand  boundary,  as  seen  in  figure  2.  The  background  neutral 
density,  n^,  is  from  some  mesurements  made  by  Pal  and  Hammer^. 
However,  this  density  is  from  a  magnetically  insulated  diode,  and  not 
a  reflex  switch.  We  use  this  value  of  nA  because,  as  mentioned 
previously,  there  is  no  hard  data  on  the  physical  parameters  inside  a 
reflex  switch.  In  additon,  see  reference  (1)  for  another  estimation  of 
the  background  neutral  density.  Now  refering  to  figure  3,  we  see  the 
code  generated  ion  flux  as  a  function  of  x,  at  time  t  =  4.54  ns  (400  x 
1.136  x  10-11  sec).  This  time  is  roughly  two  orders  of  magnitude 
greater  than  the  mean  collision  time  x,  estimated  by 


X  =  X  /  v 


v  =  ^  <E.>/m 

l 


x=x 


max 


where  <E^>  is  given  by  equation  (14).  This  means  that  we  can  assume 

that  we  have  reached  an  equilibrium  state.  The  reason  the  ion  flux  is 

depressed  near  x  =  0  is  due  to  the  fact  that  we  inject  and  push 

particles  before  data  is  collected  and  plotted.  The  dashed  line  is 

the  constant  ion  flux  F,  to  be  used  in  the  analytical  results, 

equations  (11),  (9),  and  (15).  Figure  4  is  the  code  generated  neutral 

flux  profile  and  the  dashed  line  is  equation  (11),  with  F  =  4.2  x  10^ 
2  -2 

/  cm  sec  and  X  =  5  x  10  cm.  Note  that  we  have  very  good  agreement 

between  theory  ani  simulation.  Figure  5  is  the  spatially  integrated 

distribution  of  neutrals,  as  a  function  of  energy  in  Kev,  in  the 

region  0  <  x  <  x  .  This  distribution  we  have  defined  as  G  .  (w) 

6  max  xint 

and  is  given  by 
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-  max 
s  j  g(w,x)dx  = 


G  .  ,(v) 
xint 


mF 

,  x  , 

(  max  ) 

f  i  w  ) 

.  Xmax  f  1  v  _  ) 

qE 

l  X  J 

l  qEx  J 

H  max 

1  2A  l  1  qEx  j 

L  M  max  j 

(16) 


Note  that  once  again  there  is  good  agreement  between  the  two  types  of 

results.  Figure  6  is  the  distribution  function  g  evaluated  at  x 

=xmax’  see  ecluation  (9).  Again  the  horizontal  axis  is  the  neutral 

energy  in  Kev.  The  analytical  results  agree  well  with  the  simulation, 

but  there  is  a  discrepency  around  100  Kev,  the  maximum  energy  a 

neutral  could  have.  One  possible  explanation  is  that  such  high  energy 

neutrals  are  rare,  so  the  statistics  are  low.  Another  probable 

situation  is  that  due  to  the  plotting  logic,  the  code  is  missing  these 

events  that  occur  near  the  edge  of  the  simulation  space.  However, 

figures  3-6  show  that  at  least  for  this  standard  benchmark  case, 

simulation  and  theory  agree  rather  well.  In  addition,  we  note  that 

NEUTRAL  gives  an  average  neutral  energy  of  25.4  Kev  at  x  =  xmax>  while 

from  equation  (15)  we  find  <E  >  =  28.4  Kev.  This  average  energy  is 

found  in  NEUTRAL  by  numerical  integration  of  the  neutral  distribution 

function  at  x=x  ,  see  equations  (13)  and  (15). 
max 


To  test  the  NEUTRAL  code  further  we  have  run  four  additional 
cases.  In  the  first  two  cases  we  changed  the  magnitude  of  the 
charge-exchange  cross  section,  while  keeping  the  rest  of  the  input 
parameters  the  same  as  in  the  reference  case,  see  above.  While  in  the 
second  case  we  varied  the  electric  field  E,  again  keeping  the  rest  of 
the  input  parameters  constant.  Figures  7-10  have  the  same  meaning  as 
before  except  a  is  now  twice  that  of  the  standard  case  i.e.,  a  =  4.0  x 
10'15  cm^.  From  figure  7  we  see  that  F  =  4.0  x  10^  /  cm^  sec.  This 
is  slightly  smaller  than  the  flux  used  in  the  standard  reference  case, 
we  will  discuss  this  point  later  on.  However,  we  notice  that  again  we 
have  good  agreement  between  the  analytic  solutions  and  the  simulation. 
More  importantly  these  results  all  seem  to  be  self  consistent. 
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Figures  11-14  are  for  the  case  where  the  cross  section  is  one  half  the 
value  of  the  benchmark  case,  meaning  c  =  1.0  x  10  cm  .  A  quick 
calculation  shows  that  the  mean  free  path  in  this  case  is  exactly 
equal  to  xmax>  the  width  of  the  charge-exchange  region.  This  is 
reflected  by  the  fluctuations  seen  in  figure  14,  but  even  so  the 
camparison  between  theory  and  simulation  is  still  good.  Also  note 
from  figure  13  that  because  Gxint(w)  Is  the  spatial  integral  of  the 
neutral  distribution  function,  it  is  much  smoother  than  g(w,x=x  ). 
In  figures  15-18,  we  set  the  applied  electric  field  E  to  one  half  the 
standard  field  value,  while  once  again  all  the  other  input  parameters 
are  the  same  as  in  the  reference  case.  From  figures  17  and  18  we  see 
that  NEUTRAL  underestimates  the  energy  distributions  near  zero  energy, 
by  roughly  10  percent  if  one  considers  the  numerical  values.  And 
since  the  applied  field  is  half  its  original  value,  we  of  course 
expect  that  the  maximum  neutral  energy  to  be  half  of  100  Kev,  or  50 
Kev.  This  is  verified  in  figures  17  and  18.  Figures  19-22  represent 
the  case  where  we  have  increased  the  electic  field  strength  to  twice 
the  standard  value.  From  figure  19  we  see  that  the  ion  flux  is 
reduced  by  a  small  fraction  as  compared  to  figure  15.  See  also  figure 
7  for  an  addtional  comparison  in  this  slight  flux  decrease.  Again  the 
simulation  results  appear  quite  good.  One  must  be  aware  that  in 
general  when  one  changes  some  input  parameter,  such  as  E,  then  one 
might  have  to  modify  other  values,  such  as  the  time  step  At,  in  order 
to  obtain  reasonalble  results.  It  seems  though,  from  figures  7-22, 
that  the  code  is  stable  against  changes  in  a  and  E,  as  long  as  these 
changes  are  not  too  large. 


It  seems,  from  the  discussion  of  the  simulation  results  above, 
that  the  NEUTRAL  code  finds  its  own  value  for  the  ion  flux  F,  and  for 
the  most  part  this  is  true.  But  since  we  are  supposedly  in  a  steady 
state  situation,  the  equilibrium  ion  flux  in  the  charge-exchange 
region  (0  <  x  <  xmax)  should  be  equal  to  the  ion  injection  rate,  given 
by 


F  = 


N.  . 
ini 

AAt 


(17) 
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The  quantity  A  is  an  area  weighting  factor  that  is  used  in  NEUTRAL  to 

determine  certain  physical  values,  such  as  the  charge  density.  One 

can  see  that  such  a  factor  is  needed  since  our  problem  is  one 

dimensional,  but  the  flux  and  density  are  2  and  3  dimensional  objects, 

respectively.  For  all  the  cases  mentioned  in  this  study,  A  was  set  to 
-2  2. 

10  cm  .  Using  and  At  as  given  in  the  reference  case  and 

equation  (16)  above,  we  find  that  theoretically,  F  should  have  the 
15  2 

value  4.4  x  10  /cm  sec.  From  the  figures  previously  given  we  see 
that  the  code  predicts  ion  fluxes  that  are  at  most  10  percent  off  from 
this  theoretical  value.  Ve  should  also  mention,  which  may  be  apparent 
at  this  point,  that  the  constant  ion  flux  read  from  the  plots  is 
somewhat  arbitrary,  but  at  the  same  time  one  can  make  a  reasonable 
estimation  of  the  steady  state  flux.  The  dependence  of  the  ion  flux  on 
the  variables  N^  and  At,  as  given  in  equation  (17),  is  shown  in 
figures  23  and  27.  In  figure  23  the  number  of  injected  particles  per 
time  step  was  increased  to  700,  this  is  a  factor  of  1.4  times  the 
standard  case  of  N^nj  =  500.  In  figure  27  we  have  increased  the  time 
step  by  a  factor  of  two,  this  should  reduce  the  flux  by  the  same 
factor  according  to  equation  (17),  and  this  is  essentially  the  case. 
Notice  that  this  same  figure  shows  that  if  one  increases  the  time  step 
without  adjusting  other  relevant  parameters,  that  larger  than  normal 
fluctuations  may  occur  in  the  data.  The  figures  24-26  and  28-30  are 
added  here  for  completness  only.  Like  before,  figures  24-26  are  the 
comparison  of  code  results  to  analytic  solutions  for  the  ion  flux 
given  in  figure  23,  and  the  same  goes  for  figures  28-30  except  they 
are  related  to  figure  27.  From  equation  (17)  we  see  that  the  ion  flux 
is  independent  of  the  cross  section  or  the  electric  field,  or  any 
other  physical  characteristic  that  describes  the  charge-exchange 
region.  As  was  mentioned  awhile  back,  we  noticed  a  slight  reduction 
in  the  steady  state  ion  flux  when  the  cross  section  or  the  electric 
field  were  increased  by  a  factor  of  2.  This  may  or  may  not  be 
significant.  In  both  cases  a  possible  general  argument  is  that  in  a 
given  spatial  cell  ions  are  being  lost  at  a  faster  rate  than  they  are 
injected,  reducing  the  overall  numerical  flux,  but  this  is  just  at 


first  glance.  However,  the  ion  flux  does  not  seem  to  vary  as  v  , 
2  avg’ 

defined  as  v  =  2EX/m,  as  one  may  initially  have  guessed, 
avg 
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Before  we  continue  we  should  mention  that  all  of  the  cases 
considered  so  far,  and  later  on  as  well,  the  ions  are  injected 
monoenergeticaly  from  the  left  hand  boundary  of  the  simulation  region, 
with  zero  velocity.  One  can  also  chose  to  inject  ions  with  a 
maxwellian  distribution  using  NEUTRAL,  if  the  need  arises.  In 
position  space,  the  particles  are  uniformly  distributed  in  the  first 
spatial  cell  when  they  are  injected  into  the  system. 

To  see  if  injecting  the  ions  with  a  small  velocity  would  change  our 
results,  we  used  the  reference  case  input  parameters,  but  now  each 
particle  was  given  the  initial  velocity  of  1.384  x  10^  cm/sec,  which 
corresponds  to  a  proton  at  roughly  1.0  ev  in  temperature.  The  results 
are  shown  in  figures  31-34,  and  should  be  compared  with  the  standard 
case,  figures  3-6.  For  the  most  part  one  can  see  there  is  no  drastic 
changes  between  the  two  runs.  Finally  we  will  consider  the  case  where 
we  decreased  the  size  of  the  spatial  mesh  cell  by  a  factor  of  0.5  from 
the  reference  standard.  As  one  can  see  from  figures  35,  making  the 
cells  smaller  produces  large  fluctuations  in  the  ion  flux  as  one  would 
expect,  but  there  are  only  small  ripples  in  the  neutral  flux  profile, 
see  figure  36.  In  figures  37  and  38  things  are  relatively  quiet,  but 
this  should  be  the  case  since  these  plots  are  independent  of  the 
spatial  mesh  cell  Ax.  For  comparison,  in  figures  39-42,  we  doubled 
the  size  of  Ax,  and  as  expected  the  ion  flux  shows  alot  less 
fluctuations.  We  note  that  the  NEUTRAL  code  uses  only  uniform  mesh 
cells  throughout  the  simulation  region.  This  concludes  our  evaluation 
of  the  charge-exchange  simulation  code  NEUTRAL. 


Summary  and  Conclusions 

Ue  have  found  exact  analytic  solutions,  assuming  a  constant 
charge-exchange  cross  section,  to  the  one  dimensional  Boltzmann 
equation  as  it  relates  to  the  charge-exchange  model  considered  by 
Prono  et  al.  These  results  were  used  to  benchmark  the  Monte  Carlo 
code  named  NEUTRAL,  which  was  written  to  simulate  ions  charge- 
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exchanging  in  a  background  gas.  This  exchange  process,  with  the 
resulting  flux  of  fast  neutrals,  is  thought  to  be  the  cause  of 
premature  gap  closure  in  a  reflex  switch.  In  benchmarking  the 
NEUTRAL  code  we  varied  one  input  parameter  while  keeping  the  others 
constant.  The  quantities  varied  were,  the  charge-exchange  cross 
section,  the  electric  field,  the  number  of  particles  injected,  the 
time  step,  the  speed  of  the  injected  particles,  and  the  spatial  mesh 
size.  In  most  cases  we  found  there  was  very  good  quantitative 
agreement  between  simulation  and  theory.  Thus  we  are  confident,  at 
least  for  the  constant  electric  field  case,  that  NEUTRAL  is  operating 
correctly. 


Future  Work 

We  intend  to  use  NEUTRAL  to  model  problems  that  are  more  physically 
relevant  to  the  charge-exchange  process  as  it  may  occur  in  the  reflex 
switch.  Such  as  a  cross  section  that  varies  with  ion  energy,  a 
background  gas  density  that  changes  with  distance,  and  an  electric 
field  that  varies  in  time.  Ue  will  also  use  the  results  of  the 
NEUTRAL  code  as  input  to  a  model,  yet  to  be  determined,  of  the  neutral 
shorting  mechanism  inside  the  vacuum  gap  of  a  reflex  switch. 
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Appendix  A 


The  solution  of  equation  (5),  is  not  as  straight  forward  as  one 
might  think,  so  we  will  present  the  derivation  here.  Vriting  down 
equation  (5)  for  convenience,  we  start  off  with, 


3f 

9w 


+ 


1  , 

^E(l>  * 


nA«)I 


6(w  -  s)mF 
qE 


1  +  f  dw 
Jo 


.V 


'l'  ' 


mF 


(Al) 


Creedon^  has  solved  this  equation  by  two  different  approaches, 
however  we  will  present  the  derivation  given  more  recently  by  Chen^. 
Multiplying  equation  (Al)  by  the  integration  factor 


e 


nAff' )dw' 


and  then  integrating  with  respect  to  w,  from  0  to  w,  we  find 


(A2) 


For  now  we  assume  that  l(w=0)  =  0,  and  we  will  prove  this  later  on. 

In  equation  (3)  of  the  text,  we  apply  the  Laplace  transform  and  find 


mF 

P 


(A3) 


where  we  have  again  changed  the  integration  variable  to  the  energy 
variable  w,  see  equation  (10)  .  From  equations  (A2)  and  (A3)  we  find 
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mF 

P 


and  since  the  quantity  in  the  square  brackets  is  a  constant  with 
respect  to  the  w  integration  we  can  immediately  write 


(A4) 


Using  this  result  in  equation  (A3)  we  find  the  final  form  for  the 
distribution  function,  namely 


Kv) 


(A6) 


Letting  e  -*  0,  we  see  that  we  obtain  equation  (7).  Notice  that  the 
small  factor  e  has  been  used  to  avoid  any  problems  at  w  *  0.  As  was 
mentioned  earlier,  we  assumed  that  f(w=0)  =  0,  and  we  will  now  show 
that  this  is  the  case.  Integrating  equation  (Al)  from  w  =  0  to  w  =  <*> 
and  canceling  like  terms  we  have 


00 

I(w=«)  -  f(w=o)  +  J  J(w)dw  =  jj!  (A6) 

Because  a  distribution  function  that  increased  as  w  goes  to  infinity 
is  unphysical,  we  have  the  usual  requirement  that  f  go  to  zero  as  w  •* 
<*>.  Now  using  equation  (A3)  in  (A6).  along  with  the  condition  at  w  = 
®,  we  see  that  l(w=0)  =  0.  Further  details  on  the  actual  derivation 
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of  equation  (Al),  as  was  mentioned  in  the  text,  can  be  found  in 
reference  (1). 
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Figure  14 
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Figure  32 
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Figure  35 
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